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Abstract 

The hydrokinetic model and the scheme of dynamical freeze out proposed earlier for 
the RHIC energy are extended to lower colliding energy for non-vanishing baryon 
chemical potential. In this case a new two-phase equation of state is applied where 
baryon-rich hadronic matter is described within modified relativistic mean-field the- 
ory. As an example, the approach is employed for analyzing pion and kaon spectra 
and interferometry radii of Pb+Pb collisions at laboratory energies 40 and 158 
AGeV. A good agreement is observed for the transverse mass dependence of the 
longitudinal radius for both energies studied. The transverse radii Rout,Rside re- 
produce experiment at Eiab =158 AGeV but at Eiab =40 AGeV there is a marked 
(20 %) difference from experimental points. This discrepancy is more noticeable in 
analysis of the ratio Rout/ Rside- Nevertheless, considered examples of application 
of the hydrokinetic model with equation of state including the first order phase 
transition are quite successful and further development of the model is required. 



1 Introduction 



An analysis of the experimental data in nucleus-nucleus collisions obtained at 
Relativistic Heavy Ion Collider(RHIC) [T] and very recently at Large Hadron 
Collider (LHC) [2] denotes hydrodynamics as an adequate tool to describe a 
behavior of the bulk of the matter created at these collisions. The reason could 



Preprint submitted to Elsevier 



13 December 2010 



be that the collective behavior of the quark-gluon plasma at the reached tem- 
peratures is similar to what an almost perfect fluid has. The hydrodynamic 
models use then the thermodynamic equation of state (EoS) reconciled with 
the Lattice QCD. The successful description of RHIC data within hydrody- 
namic approaches with artificial sharp freeze-out [3] as well as with dynamical 
freeze-out in the HydroKinetic Model (HKM) [^f5|6] were done and the predic- 
tions for the spectra and femtoscales at LHC are in quite reasonable agreement 

As for the applicability of hydrodynamic approach to lower energies, such as 
the SPS ones and planning experiments NICA and FAIR, it is a questionable 
issue. First, as expected, created fireballs spend only small part of its lifetime 
if any in quark-gluon phase. An evolution of the hadronic matter could be non- 
equilibrated, neither thermally, nor chemically. Since in the HKM model |1] 
the evolution and gradual decay of the hadronic matter are described as non- 
equilibrated process, it seems reasonable to try to apply this model to relatively 
small collision energies. It is the main aim of our work. 

One of the crucial points is the equation of state at large baryonic chemical 
potentials, which are typical for these energies. As known, because of the prin- 
cipal technical problem the EoS cannot be found from the lattice QCD calcula- 
tions at large baryonic chemical potentials. Currently such an estimate should 
be obtained rather from a detailed study of strongly interacting hadronic mat- 
ter supplemented by the model of quark-gluon phase at large fj,B and fj,s- In 
this paper we will use the EoS taking into account a possible phase transi- 
tion between quarks-gluons and hadrons including in-medium modification of 
hadron properties in hot and dense nuclear surrounding [5] 



2 Basic features of the model 

Let us briefly describe the main features of the HKM [5f6] . It incorporates hy- 
drodynamical expansion of the systems formed in A+A collisions and their dy- 
namical decoupling described by escape probabilities. The basic hydro-kinetic 
code includes decays of resonances into expanding hadronic chemically non- 
equilibrated system and, based on the resulting composition of the hadron- 
resonance gas at each space-time point, provides the equation of state (EoS) 
in a vicinity of this point. The obtained local EoS allows one to determine the 
further evolution of the considered fluid elements. The complete picture of the 
physical processes in central Pb + Pb collisions encoded in calculations is the 
following. 
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2.1 Initial conditions 



Our results are all related to the central rapidity slice where we use the boost- 
invariant Bjorken-like initial condition. We consider the proper time of ther- 
malization of quark-gluon matter as the minimal one discussed in the lit- 
erature, To = 1 fm/c. The initial energy density in the transverse plane is 
supposed to be Glauber-like [9], i.e. is proportional to the participant nucleon 
density for Pb-|-Pb (SPS) and Au-|-Au (RHIC, LHC) collisions with zero im- 
pact parameter. The height of the distribution - the maximal initial energy 
density e(r = 0) = eo is a fitting parameter. Also, at time r,, there is peripheral 
region with relatively small initial energy densities: e(r) < 0.5 GeV/fm^. This 
part of the matter (" corona" ) have no chance to be involved in thermalization 
process [10]. Thus, it should be considered separately from the thermal bulk 
of the matter and should not include it in hydrodynamic evolution. By itself 
the corona gives no essential contribution to the hadron spectra for central 
collisions [TU] . 

At the time of thermalization, tq = 1 fm/c, the system has already developed 
collective transverse velocities [TTP^ . The initial transverse rapidity profile is 
supposed to be linear in transverse radius r-p: 



Here a is a second fitting parameter. Note that the fitting parameter a should 
also include a positive correction for underestimated resulting transverse flow 
since in this work we do not take directly into account viscosity effects [13] 
neither at the QGP stage nor at hadronic one. In HKM formalism [5] viscos- 
ity effects at the hadronic stage are incorporated in mechanisms of the back 
reaction of particle emission in hydrodynamic evolution which we ignore in 
current calculations. Since the corrections to transverse flows which depend 
on unknown viscosity effects are unknown, we use a fitting parameter a to 
describe the "additional unknown portion" of flows, caused by three factors: 
development of the pre-thermal flows, viscosity effects in quark-gluon plasma 
and in addition event-by-event fluctuations of the initial conditions. The lat- 
ter also leads to an increase of the "effective" transverse flows, obtained by 
averaging at the final stage, as compared to the results based on the initial 
conditions averaged over initial fluctuations [H]. Since we use the last type of 
initial conditions, it should also lead to an increase of the effective parameter 
a. 




(1) 
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2. 2 Equation of state 



We match the equation of state at the point of the chemical freeze out T^^ 
which separates chemically equilibrated and non-equilibrated thermodynamic 
zones. As to the chemically equilibrated region T > T^h the EoS is constructed 
on the basis of the modified relativistic mean-field model with Scaled Hadron 
Masses and Couplings (SHMC) piTS] . In this SHMC model the Lagrangian 
density of hadronic matter is a sum of several terms: 

= -Cbar + -CmF + -Cex ■ (2) 

The Lagrangian density of the baryon component interacting via a, u mean 
fields is as follows: 



6e{bar} 



(3) 



The considered baryon set {b} consists of all baryons and low-lying resonances 
with the mass rub <1700 MeV, including strange and double-strange hyperons 
as well as antiparticles. The used cr-field dependent effective masses of baryons 
are [81I15II16I 



ml/rrib = ^bix^cr) = 1 - g^b Xa cr /nib , be {b} . (4) 

In Eqs. ([3]), (jlj) g^b and g^b are coupling constants and Xo-(o")) Xw(cr) are cou- 
pling scaling functions. 

The a-, w-meson mean field contribution is given by 



-MF — — 



b c 

uJi^u =df,uju - duujf, , U{xa(^) = 'm%{-f + -f^), f = gaN Xa (^/rriN . 

There exist only a and uq mean field solutions of equations of motion. The 
mass terms of the mean fields are 



'ml,/mm=\^miXaCr)\, {m} = a,uj. (6) 

The dimensionless scaling functions and as well as the coupling scaling 
functions Xm, depend on the scalar field in combination Xo-(o") cr. Following 



4 



|16] we assume approximate validity of the Brown-Rho scaling ansatz in the 
simplest form 

$ = $Ar = = $^ = = 1 - /. (7) 

The third term in the Lagrangian density ([2]) corresponds to meson quasiparti- 
cle excitations, where particles with < 1100 MeV are included The choice 
of parameters and other details of the SHMC model can be found in piTS] . 
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Fig. 1. The temperature dependence of effective masses of the nucleon, uj and p 
excitations (solid line) and of the cj-meson excitation (dashed line) calculated within 
the SHMC model for two values of the baryon density ne- 

Within SHMC model different thermodynamical quantities in thermal equilib- 
rium hadron matter are calculated at fixed temperature T and baryon chemical 
potential /xb- In Fig H] the T and /xb dependence of hadron masses is pictured. 
Note that these masses sharply decrease in the vicinity of the expected phase 
transition and thereby affect on the fireball evolution. 

The quark-gluon plasma phase is treated in the conventional quark bag model 
with massive quark and gluons. A two-phase model allowing for the first order 
phase transition from the hadron phase (H) to the strongly coupled quark 
gluon plasma (Q) is constructed by means of the Gibbs conditions. These 
are conditions for thermal {T^ = T^), mechanical {P^ = P^) and chemical 
(/^B = /^B, /^Sr = ^J^ftv) equilibrium. 

There is a reasonable agreement of our two-phase thermodynamics with avail- 
able QCD lattice calculations [81IT5] . For /x^ =0 some discrepancy is seen in 
energy density near the critical temperature Tc due to the fact that by con- 
struction we have the first order phase transition while in lattice calculations 
we deal with a sharp crossover-type transition. The colliding energies of inter- 
est are below an anticipated critical point and the first order phase transition 




5 




/ig, MeV GeV/fm' 

Fig. 2. Left panel: Isoentropic trajectories in the T — fj,B plane for two values of 
entropy-to-baryon number ratio. The quark-hadron phase boundary [15] and the 
freeze out curve [17] are shown by the dotted and dashed dash-dashed lines, re- 
spectively. Experimental points are from review article [T7j . Right panel: Equation 
of state within two-phase model for the same values of the entropy baryon number 
ratio as in the left panel. 

is expected. Unfortunately, there is no lattice data for these large /i^ range. In 
the left panel of Figj2]we demonstrate isoentropic trajectories for two fixed ra- 
tios of entropy S to baryon number A^^, S/Nb =19.2 and 32.8, corresponding 
to that of an expanding fireball at Eiab = 40 and 158 AGeV. When a trajec- 
tory reaches the phase boundary it is flattening in accordance with the mixed 
phase of the Maxwell construction in the two-phase model (see the right panel 
in Fig|2]). Assuming chemical equilibrium our hydrodynamics is used till the 
chemical freeze out curve shown in Fig|2]by the dotted line. At temperature 
below Tch we apply non-equilibrium treatment. 

The EoS for baryon-rich fireballs described above is used up to the tempera- 
ture of the chemical equilibrium Tch- We should match this chemically equi- 
librated region T > Tch and non-equilibrated hadronic region at T < Tch- 
The last is described in the HKM as an ideal gas of hadrons and their res- 
onances in correspondence with the concept of chemical freeze-out. We fix 
the chemical freeze-out temperature and baryon chemical potential to be 
{T,//b} = {160 MeV, 220 MeV} for ^Lab = 158 GeV collisions and {T,/xb} = 
{148 MeV, 370 MeV} for ^Lab = 40 GeV coUisions, based on Tch{y/s]^), 
/^bIa/s/wv) parameterizations obtained from particle number ratio analysis 
performed in [19] . In both cases, strange chemical potential is calculated from 
the condition of zero strange density. The proper choice of Tch and guar- 
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antees us the correct particle number ratios for all quasi-stable particles (here 
we calculate only pion and kaon observables) . Below Tch a composition of 
the hadron gas is changed only due to resonance decays into expanding fluid. 
We include = 326 well-established hadron states made of u, d, s-quarks 
with masses up to 2.6 GeV from PDG table [18], as well as the most possible 
decay channels of them. The EoS in this chemically non-equilibrated system 
depends now on particle number densities rij of all the Nh particle species 
i: p = p(e, {rij}). Since the energy densities in expanding system do not di- 
rectly correlate with resonance decays, all the variables in the EoS depend 
on space-time points and so an evaluation of the EoS is incorporated in the 
hydrodynamic code. We calculate the EoS below Tch in the Boltzmann approx- 
imation of ideal mult i- component hadron gas. In addition to standard HKM 
method we add now excluded volume Vq (equal for all hadrons) which may be 
important for baryon-rich systems, according to the following formula: 



the superscript id refers to ideal gas case. We solve this equation numerically, 
and obtain other thermodynamic quantities using general thermodynamic re- 
lations and the solution of Eq. . 

We match energy and baryon number densities in both EoS at chemical 
freeze-out point, and minimize pressure difference there as well. We choose 
the excluded volume value Vq = 0.2 fm^ from the condition of the best 
matching. The corresponding energy and net baryon densities are {e, ne} = 
{0.4 GeV/fm^, 0.085 l/fm^} for ^Lab = 158 GeV and {e, hb} = {0.34 GeV/fm^, 
0.12 1/fm^} for -^Lab = 40 GeV collisions, which correspond to the chemical 
freeze-out parameters shown above. The chosen excluded volume in chemically 
non-equilibrated ideal hadron gas thus mimics more complicated interactions 
happening in hadron matter in two- phase model before chemical freeze-out. 

2. 3 Evolution 

At the temperatures higher than Tch the hydrodynamic evolution is related 
to the possible quark-gluon, mixed and hadron phases which are in chemical 
equilibrium with some baryonic chemical potential. The evolution is described 
by the conservation law for the energy-momentum tensor of perfect fluid: 



Since in perfect hydrodynamics the entropy is conserved, the system evolves 
along trajectories shown in Figj2l 



p(T, {/i,}) = ^pf(T, {fi^-vo-p}) 





(9) 
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At T < Tch the system evolves as chemically non-equilibrated hadronic gas. 
The concept of the chemical freeze-out implies that afterwards only elastic 
collisions and resonance decays take place because of relatively small densi- 
ties allied with a fast rate of expansion at the last stage. Thus, in addition 
to Eq.Q, the equations accounting for the particle number conservation and 
resonance decays are added. If one neglects the thermal motion of heavy res- 
onances, the equations for particle densities ni{x) take the form: 

df,{ni{x)u''{x)) = -TiUi^x) + ^ bijTjnj{x) (10) 

j 

where bij = BijNij denote the average number of i-th particles coming from 
arbitrary decay of j-th resonance, Bij = Tij/Tj^tot is branching ratio, Nij is a 
number of i-th particles produced in j — )■ z decay channel. Eqs.dH]) and Eqs. fllOl) 
are solving simultaneously with calculation of the EoS, p{x) = p(e(x), {ni{x)}) 
at each point x. 



2.4 System's decoupling and spectra formation 

During the matter evolution, in fact, at T < Tch, hadrons continuously leave 
the system. Such a process is described by means of the emission function 
S(x,p) which is expressed through the gain term for i-th sort of particles, 
Gi{x,p), in Boltzmann equations and the escape probabilities [HIS] 

Si{x,p) = Gi{x,p)Vi{x,p) 
00 

Vi{x,p) = exp{- J dsRi+h{s,r + ^{s - t),p)) . 

For particle emission in the relaxation time approximation Gt^ ~ f-KR-n+h + 
Gh^tt where f-n^x^p) is the Bose-Einstein phase-space distribution, Ri^h{x,p) 
is the total collision rate of the particle, carrying momentum p, with all the 
hadrons h in the system in a vicinity of point x and the term Gn^i describes 
an inflow of particles into phase-space point {x,p) due to the resonance de- 
cays. It is calculated according to decay kinematics with simplification that 
the spectral function of the resonance H is — {mH)"^). The cross-sections 
of interacting particles in the hadronic gas, which determine the escape prob- 
abilities V{x,p) and emission function S{x,p) via the coUision rate Ri+h, are 
calculated in accordance with the UrQMD method ^U\. The spectra and cor- 
relation functions are found from the emission function 5* in the standard way 
(see, e.g.,Ilj). 
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3 Results and conclusions 



The results for transverse mass spectra and interferometry radii are presented 
in Fig. [3] for the two SPS energies. The two fitting parameters: maximal initial 
energy density eo and the initial transverse velocity (associated with parameter 
a) averaged over the energy profile (vt) are equal to eo = 5 GeV/fm^ ((e) = 
3.5 GeV/fm^) for ELab=^0 AGeV and eo = 8 GeV/fm^ ((e) = 5.6 GeV/fm^) for 
-ELab=158 AGeV and {vt) = 0.196 for both energies. As one can see the pion 
and kaon transverse mass spectra, their slopes as well as the absolute values 




Fig. 3. Pion and kaon spectra (top), pion Rout, ^sidc (middle), -Riong and -Rout/-Rside 
ratio (bottom) as a function of transverse momentum at the given centrality. Exper- 
imental points for transverse spectra and interferometry radii are taken from |21j 
and |22|23] . respectively. 
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are in a quite reasonable agreement with available experimental data |21j . 



A good agreement with experiments [22]l23j is also observed for the longitudi- 
nal radius Riong for the both energies studied. The Rout and Rside interferom- 
etry radii are well reproduced for the top SPS energy including even the ratio 
Rout/ Rside which is more sensitive to model parameters. For Eiab =40 AGeV 
there are some deviations of the calculated curves from experimental points 
which are more pronounced for the Rout/ Rside ratio. 

Summarizing, one can note that the hydrokinetic model can serve as the basic 
dynamical model also for low-energy collisions where the equation of state 
corresponding to the first order phase transition should be used rather than 
that with crossover transformation between quark-gluon and hadron matter. 
Such an EoS at large chemical potentials can be successfully constructed by 
matching the relativistic mean-field approximation for EoS at T > T^h with 
ideal hadronic gas with hadron masses below 2.6 GeV including the excluded 
volume assumption at T < Tch- This approach reasonably reproduces not only 
bulk properties of produced particles but also such delicate characteristics as 
interferometry radii. It is worthy to note that even at Eiab= 40 AGeV the 
description of experimental data requires to take into account quark-gluon 
degrees of freedom in initial thermodynamic conditions. The HKM may be 
useful in predictions of collective observable in the energy considered in future 
FAIR and NICA experiments. 
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